Improper ferroelectricity and multiferroism in 2H-BaMn03 
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Using first-principles calculations, we study theoretically the stable 2H hexagonal structure of 
BaMnOa. We show that from the stable high temperature PGa/mmc structure, the compound 
should exhibit an improper ferroelectric structural phase transition to a PGscm ground state. Com- 
bined with its antiferromagnetic properties, 2H-BaMn03 is therefore expected to be multiferroic at 
low temperature. The phase transition mechanism in BaMnOg appears similar to what was reported 
in YMnOa in spite of totally different atomic arrangement, cation sizes and Mn valence state. 

PACS numbers: 63.20.kk, 78.30.-j, 63.20.dk, 75.85. -|-t, 63.20.-e 



Multiferroic compounds combining ferroelectric and 
(anti-)ferromagnetic orders present both fundamental 
and technological interests and have been the topic of 
intensive researches during the last decadei"— . The iden- 
tification of new single-phase multiferroics has been par- 
ticularly challenging. In this search, a special emphasis 
was placed on the family of multifunctional ABO3 com- 
pounds within which belong many popular ferroelectric 
and magnetic materials. Unfortunately, few multiferroic 
ABO3 oxides have been identified to date. This scarcity 
was justified in terms of typical antagonist metal d-state 
occupacy requirements^: on the one hand, ferroelectric- 
ity in typical compounds like BaTiOs arises from the hy- 
bridization between occupied O 2p orbitals and empty Ti 
3d states and requires d^-ness while, on the other hand, 
magnetism requires partial d orbitals occupacy, so that 
ferroelectricity and magnetism often appear to be mutu- 
ally exclusive. 

Exceptions to this rule exist however. In the antiferro- 
magnetic CaMnOa perovskite, a ferroelelectric instability 
exists in spite of partial d-state occupacy. In bulk sam- 
ples, this instability is suppressed by antiferrodistortive 
oxygen motions but it was predicted theoretically^, and 
recently confirmed experimentally^, that CaMnOa can 
be made ferroelectric, and de facto multiferroic, under 
appropriate tensile epitaxial strain. In the same spirit, 
it was proposed that, due to its larger volume, the per- 
ovskite form of BaMnOa should develop a ferroelectric 
antiferromagnetic ground state^ii. 

BaMnOa does however not naturally crystallize in the 
perovskite structure with corner-sharing oxygen octahe- 
dra. Due to its large Goldschmidt tolerance factor, it 
stabilizes instead in a hexagonal 2H structure with face- 
sharing oxygen octahedra (Figure [1])^'^. Increasing pres- 
sure decreases the hexagonal character, yielding at high 
pressure various possible forms including a 411 hexagonal 
structure combining corner- and face-sharing octahedra, 
which is intermediate between 2H and cubic perovskite 
structures^. The stability of the 211 structure was con- 
firmed from first-principles calculations^^ but its exact 
symmetry remains under debate. 

From the similarity of spectra with BaNiOy^, Hardy 
originally assigned the 2H phase of BaMnOa at room 
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FIG. 1. Crystallographic structure of 2H — BaMnOa com- 
pound at low temperature {P6aC7Ji space group), (a) Face- 
sharing MnOe octahedra forming chains along the c axis, (b) 
Triangular arrangement of the MnOe octahedra in the (a, b) 
plane. 



temperature to the polar PQ-^mc space group^. Com- 
bined Raman and infra-red phonon measurements per- 
formed by Roy and Budhanii^ are compatible with such 
a polar character, that is also in line with a recent report 
of room temperature ferroelectricity in 2H — BaMnOa by 
Satapathy et al^^. 

Nevertheless, Christensen and Ollivier reported that 
2II-BaMn03 should better belong to the PQ^/mmc space 
group>i^. We notice that a similar re-assignment to 
PQ^/mmc was also proposed for BaNiOg^^ that had 
served as reference material to Hardy. This was also 
the choice of Cussen and Battlo^S, who found no clear 
evidence to assign the room temperature structure of 
2H-BaMn03 to the polar PQ^mc space group but in- 
stead chose the non-polar PQ^/mmc space group. Go- 
ing further, they highlighted a structural phase transi- 
tion between room temperature and 80 K to a polar 
P6acm space group, yielding unit-cell tripling. More- 
over, while Christensen and Ollivier— were not observ- 
ing any clear magnetic order above 2.4 K, Cussen and 
Battlei^ reported an antiferromagnetic ordering below 
Ttv = 59-ftr in the PQ^cm phase, with neighboring spins 
antiferromagnetically coupled along the polar axis and 
forming a triangular arrangement perpendicular to it. 
Although it was not explicitly emphasized, according to 
this, BaMnOa should be multiferroic at low tempera- 
ture in its stable 2H polymorph. Nevertheless this has 
to be taken with caution since Cussen and Battle can- 
not rule out the possibility of a centrosymmetric P3cl 
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space group (subgroup of PG^/mcm) at low tempera- 
turei^. Also, no consensus has been achieved yet regard- 
ing the stable structure at room tcmperaturei^. 

In order to get additional insight into the possible 
ferroelectric nature of 2H — BaMnOs, we performed 
first- principles calculations. We show that the high- 
temperature phase should belong to the P63 / mmc space 
group and that the compound should undergo an im- 
proper ferroelectric structural phase transition to a 
PG^cm phase, through a mechanism similar to what hap- 
pens in hexagonal YMnOa. Combined with its anti- 
ferromagnetic properties, this study confirms that 2H — 
BaMnOs should be multiferroic at low temperature. 

Our calculations have been performed within density 
functional theory (DFT) using the CRYSTAL09 pack- 
ageii. We used the Bl-WC hybrid functionali^ that was 
previously found appropriate for the study of multifer- 
roic oxidesi^. Rclativistic pseudopotcntials and a 3^ va- 
lence basis sets were used for the barium atoms2£. Small 
core Hay-Wadt pseudopotentials, associated with 2( va- 
lence basis sets^i, were used for the manganese atoms. 
All-electron 2^ basis sets, specifically optimized for 0^~, 
were used for the oxygen atoms2^. The calculations were 
carried out on a 4 x 4 x 8 fc-point grid for the tripled unit 
cell (30 atoms). The polarization was computed from 
the displacement of Wannier function center o^'^'^^ . We 
worked within a collinear-spin framework assuming an A- 
type antiferromagnetic order with spins antiparallcl along 
the c direction and aligned in the (a, b) plane. Calcula- 
tions with a ferromagnetic ordering provided similar re- 
sults. We performed geometry optimization of atomic po- 
sitions until the root mean square values of atomic forces 
were lower than 5 x 10~^ Ha.Bohr"^ and until the root 
mean square of atomic displacements were lower than 
5 X 10~^ Bohr. For better comparison with the experi- 
ment, results reported below have been obtained at fixed 
experimental lattice parameters as reported by Cussen 
and Battle at room tcmpcraturoi^. Similar results have 
been obtained adopting the lattice parameters reported 
at 1.7 K and 80 K. Full structural optimization including 
lattice parameters were also performed providing quali- 
tatively similar results but yielding a smaller ferroelectric 
distortion, due to small theoretical errors on the lattice 
parameters. 

First, we performed a structural optimization of 2H — 
BaMnOs starting from a polar PQ^mc configuration as 
initially proposed by Hardy^. Surprisingly, the system 
went back to a PQ^/mmc configuration, with atomic po- 
sitions (one degree of freedom, xq — 0.1473) in good 
agreement with what was reported by Cussen and Battle 
{xq = 0.1495)i^. This demonstrates theoretically that 
the system does not want a priori to be ferroelectrically 
distorted but prefers to adopt a non-polar Pd^/mmc 
symmetry as first proposed by Christensen and Ollivier. 
This is in confiict with the recent claim of a room tem- 
perature PG^mc ferroelectric phase by Satapathy et a&^. 
However the crystal structure they report corresponds in 
fact to a P63 /mmc phased and is therefore inconsistent 
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TABLE I. Atomic positions (in fractional coordinates) in the 
P63cm structure of 2H-BaMn03. Experimental values are 
those of Cussen and Battle at 1.7K and 80K— . 



Ki K3 total 



Exp. (1.7K) 0.015 0.144 0.154 0.534 0.574 
Exp. (80K) 0.035 0.139 0.035 0.420 0.445 
Theory 0.001 0.034 0.005 0.766 0.767 



TABLE IL Amplitude (in A) of symmetry allowed displace- 
ments for experimental and theoretical structures (at fixed 
lattice parameters). 



with the polarization they measure. In order to further 
confirm our result, we computed the zone-center phonon 
modes of the relaxed PQ^/mmc phase without identify- 
ing any zone-center unstable mode, the lowest polar T2 
mode being at a frequency around 56 cm^^. 

In order to test the possibility of a phase transition to 
a PQ^cm or a P3cl phase, we computed the phonons at 
the zone-boundary point q=(l/3, 1/3, 0). In this case, 
we identified an unstable mode at 25i cm~^ of A'3 sym- 
metry. This mode is related to oxygen and manganese 
motions along the polar axis. Its condensation produces 
a unit-cell tripling and brings the system into the PGacm 
symmetry. We notice that no unstable mode of Ki sym- 
metry is identified, rulling out the possibility of a transi- 
tion to a phase of PQ^/mcm or P3cl symmetry. 

Atomic relaxation of the PQ^cm phase has been per- 
formed, yielding a gain of energy of 21 meV/f.u. respect 
to the PQs/mmc structure. The relaxed atomic positions 
are reported in Table |T] and show qualitative agreement 
with experimental data. BaMnOa appears in our hybrid 
functional calculation as an insulator with an electronic 
bandgap of 4.15 eV. As allowed by the polar nature of 
the PGacm space group, this phase also develops a siz- 
able spontaneous polarization Pg = 0.92 /xC.cm"^ along 
the polar axis. Since this polarization can be either up 
or down (and is so a priori switchable), this phase can 
be labelled as ferroelectric. 

Interestingly, although the condensation of the 
mode breaks inversion symmetry, this mode is non-polar 
(i.e. its mode effective charge is zero) and its condensa- 
tion alone does not induce any significant polarization^^. 
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FIG. 2. Evolution of tlie energy (in meV) of 2H-BaMn03 in 
terms of the amplitude Qk^ and Q„- (in fractional units) of 
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the ifa and modes. The P&z/mmc structure was taken 
as reference. 



Further insight into the mechanism yielding a significant 
polarization at the phase transition from PQ^/mmc to 
Pdscm is therefore needed and can be obtained from a 
decomposition of the atomic distortion condensed at the 
transition into symmetry adapted modes, as achievable 
with the Amplimodes software of the Bilbao Crystallo- 
graphic server^^^. 

The results of this decomposition are reported in Ta- 
ble ini Both theory and experiment find the structural 
distortion to be principally dominated by a K3 char- 
acter. Then, and although underestimating the exper- 
iment, another significant atomic motion observed theo- 
retically has a r2 character. As further discussed below, 
it is in fact this additional polar distortion that will pro- 
duce a sizable polarization. We notice that the appear- 
ance of additional T^ and A'l atomic distortions are also 
possible by symmetry but arc not observed in our calcula- 
tions. This last result contrasts with experimental data 
that suggest the presence of a sizable Ki contribution. 
We notice however that experimental data at 1.7K and 
80K are not fully compatible in that respect, so that fur- 
ther refined experimental characterization of this phase 
would likely be very useful to clarify that issue. 

The interplay between dominant K3 and atomic 
distortions can be better understood from an expansion 
of the energy around the P63 / mmc phase in terms of 
their amplitude Qks and Qp- : 



primary mode. Then the sizable polarization arises from 
the additional condensation of the stable polar F^ sec- 
ondary mode which, as further illustrated in Fig. [31 re- 
sults from a progressive shift of the minimum of the single 
weU, produced by Qk^ through the /^siQ/^^Qp- coupling 
term. 
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FIG. 3. Energy (in meV) as a function of Fj allowed dis- 
placements at fixed Qk'j (in fractional units). 

This mechanism is similar to what was reported by 
Fennie and Rabe — in YMnOa and allows us to iden- 
tify 2H-BaMn03 as an improper ferroelectric. The evo- 
lution of the polarization with Qk^ is represented in 
Fig. m We first notice that the polarization esti- 
mated within the model from the amplitude of Qp- as- 
suming constant Born effective charges {Z^^ = +2.10e. 



~ +6.4:5e and Zq = — 2.85e along the caxis in the 
PQy, / mmc phased) properly reproduces the DFT calcu- 
lation. As discussed by Fennie and Rabe—, two indepen- 
dent regimes can be identified: at small Qk^ amplitude, 
Qp- (and therefore Pg) evolves like Q\,^ while at larger 
amplitude it evolves linearly. Contrary to YMnOa that 
has its ground state in the linear regime, BaMnOs is still 
in the small amplitude regime in our calculation. 

Although they both crystallize in the same hexagonal 
P63 / mmc space group at high-temperature, the similar- 
ity between BaMnOs and YMnOa is astonishing. First, 
the valence state of the cations are different in both com- 
pounds {+2/ -f 4 in BaMnOa versus +3/ + 3 in YMnOa). 
The cation sizes are also distinct in both compounds 



■F{Qk,, Qr-) = a2oQi', + "02 Qp- + /34oQk + (^oaQ 



%2 Q/f, Qp- 



(1) 



The coefficients of this energy expansion were deter- 
mined from DFT calculations on a y/?>apE x VSbpE xcpE 
optimized PQs/mmc supercell with various amplitudes of 
Q^3 and Qp- condensed. The results are summarized in 
Fig. [2] and are consistent with what was previously dis- 
cussed. The F^ polar mode is stable and associated to 
a single well while the K3 mode is unstable and has a 
typical double- well shape. 

The phase transition of 2H-BaMn03 from PQ^/mmc 
to PQ^cm is therefore clearly driven by the unstable 



DFT 

model: «(a„,/p„)" 
model: Q^«(ajfi,j" 



P-(-p,/(2a„,))'"Q^-' 
P-(-|3„/(2(ij^))"'Q^ 



U 



1 
(fractional units) 



FIG. 4. Spontaneous polarization (in fj,C.cm ^) as a function 
of the primary order parameter Qks (in fractional units). 
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yielding a Goldschmidt tolerance factor t > 1 in BaMnOs 
and i < 1 in YMnOs. Consequently, the atomic struc- 
ture is totally different in both compounds: the Mn 
atom is surrounded by face-sharing oxygen octahedra in 
BaMnOs and corner-sharing oxygen trigonal bipyramids 
in YMnOa. The atomic motion associated to the un- 
stable mode is also rather different. In YMnOa, K-^ 
motions correspond to a tilt of MnOs bipyramids accom- 
panied by a shift of Y atoms and O atoms forming the ba- 
sis of MnOs bipyramids along the polar axis. This leads 
to a more symmetric environment for the Mn atom, re- 
ducing its oxidation degree and making the charge of all 
O atoms similar. In BaMnOa, the mode is associated 
only to a shift of Mn and O atoms along the polar axis, 
increasing the absolute value of the oxidation degree of 
all atoms. The fact that both compounds exhibit an un- 
stable mode, that can drive improper ferroelectricity 
appears therefore rather fortuitous. 

Finally, we studied the magnetic properties of ground 
state 2H — BaMnOs. In their seminal work, Christensen 
and OUivier— did not identify any clear magnetic order 
above 2.4 K but more recently, Cussen and Battle^ re- 
ported an antiferromagnetic ordering below T^r = b^K in 
the Pdscm phase: they proposed a configuration compa- 
rable to YMnOs in which neighboring spins are antifer- 
romagnetically aligned along the polar axis and forming 
a triangular arrangement perpendicular to it. 

Although working within a collinear-spin approxima- 
tion, we estimated the coupling constants Jij between 
spins at site i and j in the PQ^cm phase through the fit of 
the Heisenberg effective Hamiltonian limited to nearest- 
neighbor interactions: Hueisenberg = - Yji<j Jij Si ■ Sj- 
2II-BaMn03 presents two distinct magnetic interactions: 
one along the polar axis and one perpendicular to it. 
We found a strong antiferromagnetic coupling Jc = 
—32.77 meV along the polar axis and a much weaker 
interchain antiferromagnetic coupling Jab ~ —0.14 meV 
in plane, comparable to that computed in YMnO^^. 
The estimate of Jab has certainly to be taken with more 
caution than Jc'. in YMnOa, collinear-spin calculations 



significantly underestimated the magnitude of Jab on 
a similar triangular arrangement of spins^i (0.59 meV 
instead of 2.3-3 meV experimentally). Nevertheless, a 
strong anisotropy of the magnetic interactions seems to 
be likely and can be explained from simple geometry ar- 
guments. Along the polar axis, there are face-sharing oc- 
tahedra and the super-exchange term between two neigh- 
boring Mn atoms goes through three different Mn-O-Mn 
exchange paths, increasing the overlap between Mn 3d 
and O 2p orbitals and the magnetic coupling. On the 
other hand, Mn-Mn interchain distances are longer and 
the super-exchange term in-plane has to go consecutively 
through O and Ba atoms (see Fig.[T]); the overlap between 
orbitals will be smaller and the magnetic exchange inter- 
actions weaker. So, although it is proposed that both 
share the same magnetic arrangement, the antiferromag- 
netic coupling in 2II-BaMn03 appears substantially dif- 
ferent than in YMnOs for which the strongest magnetic 
interactions are in-plane. 

In summary, we propose in this Letter that 211- 
BaMnOs is an improper ferroelectric. Our calculations 
moreover confirms its antiferromagnetic character, mak- 
ing it a multiferroic compound amazingly similar to 
YMnOa , in spite of a radically different atomic arrange- 
ment. Improper ferroelectrics have recently generated 
increasing interest in view of their unusual electrical 
properties'^. Not so many ABO3 improper ferroelectrics 
have been reported yet and our study illustrates that 
such behavior can also happen in hexagonal 2H com- 
pounds. We hope that our theoretical work will motivate 
the search of new improper ferroelectrics and multifer- 
roics in this class of compounds. 
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